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Introduction 


Explorer 35, injected into selenocentric orbit on July 22, 1967, 
experiences a decrease in its moment of inertia as it passes through the 
lunar shadow. This is apparently due to the cooling of the spacecraft 
and most probably inward contraction of the solar paddles. Because angular 
momentum is conserved, the spin period decreases slightly (a few parts per 
thousand) and the psuedo-sun pulse generated onboard the spacecraft no 
longer accurately indicates the sunward direction. This report describes 
an empirical approach to solving the resulting satellite aspect problem. 

A variable spin period of the form 

T = a + be' t/t: ° (1) 

is assumed. The measured spin periods during the recovery following a 
shadow passage show good agreement (See Figure 1) with this analytic form. 
During the shadow the absence of a true sun pulse precludes the accurate 
determination of the spin period. However, to the extent that the azimuth 
angle of the magnetic field vector in the satellite spin plane remains con- 
stant throughout the shadow, the magnetic field measurements themselves 
give an estimate of the satellite aspect. A comparison of the apparent 
rotation of the magnetic field with the apparent rotation calculated using 
the methods to be described is shown in Figure 3. This plot and others 
like it for subsequent shadow passes attest to the accuracy of the calcu- 
lations using a time dependent spin period of the form in 1 above. 

In order that the method of analysing the spacecraft aspect in the 
shadow be as clear as possible it seems helpful to present a brief des- 
cription of the pertinent facts of the spacecraft operation at this point. 
The direction of the sun relative to the spacecraft coordinate system is 
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measured once every four telemetry sequences (1 sequence = 81.808 seconds). 
This is measured ,by counting the time from a fixed reference point in 
the telemetry sequence to the next "see sun" pulse from the onboard sun 
sensor (for a discussion of the telemetry format and optical aspect 
system see AIMP (IMP-D) Technical Summary Description, 1967). This time 
(referred to as the "sun time") together with the time between consecutive 
sun pulses (the raw spin period) are telemetered once every four sequences. 
This enables one to calculate accurately the azimuthal orientation of 
the spacecraft about the spin axis at any time. 

Such an accurate determination of the spacecraft azimuth requires 
the calculation of a refined spin period. This is done by calculating 
the time between successive sun times and dividing by the integral number 
of revolutions made by the spacecraft. The number of . revolutions , N, 
is calculated using the raw spin period. Thus the refined spin period 
is given by 

T - c ° + (2) 

N 


where C 0 = the time between successive sun sightings. 

C^="sun time" at beginning of interval. 

C 2 = "sun time" at end of interval. 

N = integral number of revolutions between sun sightings. 
N in turn is given by 


N + A = 


C 0 + C 2 - C 1 


S 


(3) 
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where S is the raw spin period and A is a fraction. The fraction A will 
be less than 0.5 as long as the time C Q corresponds to less than about 
40 sequences. Thus N is easily found by rounding the quantity (N + A) 
to the nearest integer. 

If the spin period remained constant through the shadow it would 
be a simple matter to accurately determine the azimuthal orientation 
of the spacecraft. Actually the spacecraft spins up in the shadow, and 
by the end of the shadow a large error (of the order of 360°) has accum- 
ulated in the satellite orientation computed from the optical aspect 
data. It is this error that must be corrected. 

The time the spacecraft enters the shadow is measured to the nearest 
half sequence except every fourth sequence when it is accurate only to 
the nearest sequence. This information is telemetered in a sun flag 
which indicates whether the sun sensor is seeing the sun or not. 

Thus the situation is as follows: The spin period and sun direction 
are known accurately at the time of the last sun sighting before the shadow. 
Then there follow 0 to 4 sequences of measurements before the satellite 
enters the shadow and the spin period begins to change. Upon entry into 
the shadow the spacecraft begins to spin up and the spin period begins to 
decrease. This decrease continues until the spacecraft again enters the 
sunlight where upon it begins to de-spin again and the spin period begins 
to increase, finally approaching its pre-shadow value. After leaving the 
shadow there will be 0-4 sequences before true solar aspect data are 
obtained. The total fractional change in the spin period is about 0.0015 
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(0.157o) over the period of the shadow, typically 40 min. (35 sequences). 

With this understanding of the problem the actual calculation of 
the three free parameters in 1 is relatively straightforward. In carrying out 
the calculations the following parameters are used. Inspection of Figure 
2 illustrates their significance. 

A tf = time between the last true sun pulse and the beginning of 
the shadow 

A t s = duration of the shadow 

A tf = time between shadow end and next true sun time 
T f = spin period before shadow 

= first spin period measured after shadow 
T 2 = second spin period measured after shadow 
Cf = last M sun time" before shadow 
= first n sun time" after shadow 
S = raw spin period 

In addition the following parameters will be used: 

Tf = spin period at the end of the shadow 
R(t) = Actual number of revolutions of the spacecraft since last 
true sun before shadow (time 0 on Figure 2) 

R 0 = R (t = Atf 4* At s + Atf) = Total number of spins made between 
true sun pulses 

N Q = Integral number of rotations of spacecraft between true sun 
pulses . 

Convenient units are as follows: 

At's and t*s in seconds 

Cf, Cf, S in 800 cps clock counts 

R, R OJ N in revolutions. 
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It is often useful to express the At f s in units of sequence length = 
81.808 seconds and t G will also be measured in units of sequence length. 
Analysis 

The problem now becomes one of fitting a spin period of the form 1 
to the satellite data. There are 3 free parameters in 1 so three inde- 
pendent equations are required to determine these parameters. Two are 
obvious: the spin period before and after the shadow must be equal to 

that measured. The third condition is that the integrated number of 
revolutions between the true sun times be equal to that achieved by the 
spacecraft. To this end Tf and R Q must be determined since they are 
not directly measured. 

First Tf may be found by using a linear interpolation from T and 
t 2 at the end of the shadow. (Earlier work also fit these data to an ex- 
ponential to find Tf, but it was found to be unnecessary, the linear 
approximation being of more than adequate accuracy for the short interval 
involved.) 

Thus 

T f = ) ^ fc f + 2°)* (2) 

A little thought will reveal that R Q is given by 

R 0 = N 0 - l f — . C .i- (3) 

s 

where N 0 is to be determined. It may be found by integration if the spin 
period as a function of time is known. In this work the spin period is 
assumed to be given by 

r ±9 t ^ At£ 

a+be - ^ -At ^/t 0 ],Ati < t £ 1^ 

Tf ^ — — ) (t-Atg-At^) ,L^ < t ^ L 2 


T (t) = 


(4) 
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where = Atf + At s and L 2 = Lf + Atf. 
Now R 0 may be obtained from: 


-J L2 ^ 


T(t) 


= Atf + At £ 


At c 


- t c In 


/a+be 


°> + 


Atf 


a +b ' T f + ( t 2-t x ) At f 


(5) 


Since initially the values of a, b, and t Q are unknown we still do 
not know R Q . However , the present goal is to obtain N Q so it is sufficient 
only to know R Q to better than + 0.5 revolutions. It is found that the 
following approximate values of a, b and t Q give R 0 to + 0.1 revolutions. 


a T£ - .0001 seconds 
b = T-j_ - a 

~ “ t^) + .0001 seconds 

t G = 24 sequences ^ 1963 seconds 

Thus by rounding R Q to the nearest integer N Q is determined. This 
is exact if the estimate of R Q was within + .5 revolutions and then 3 
may be used to determine R Q accurately. (The accuracy of N Q is easily 
verified since an error shows up in the magnetic field data as an apparent 
rotation of one or more complete revolutions.) With the information derived 
above it is a relatively straightforward matter to set up the three con- 
ditions which determine the three unknown parameters in 1, a, b and t . 

The first two equations are obtained by requiring that the spin 
period be continuous at the beginning and end of the shadow. When the 
satellite enters the shadow t = At^ and T = so that 


a + b = 


( 6 ) 
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At the end of the shadow t = and t = Tfo Thus 

a -f be “( At s/ t o) = Xf . 


( 7 ) 


The third and final condition is set by requiring that the predicted 
number of revolutions between the last real sun pulse before the shadow 
and the first real sun pulse after the shadow is equal to the actual 
number of revolutions, R 0 . This expressed by 5* Since some of the terms 
in 5 involve only the known parameters (i e e 0 , do not involve a, b, or t Q ) 
it is useful to define another quantity, R s , the number of revolutions 
the spacecraft makes during the actual shadow » 


R s = 


R„ 



At^ 

T f Cl i h ) At 

4 r 


(9) 


Then the third condition may be rewritten as: 

At ► , . -(At s /t 0 ). 

At s _ £o in + be ) 


Rs « 


a 4- b 


(8a) 


Since the change in the spin period during the shadow is observed 

to be only about one one -thousandth of the spin period it is easily seen 

that -g ~ .001 « 1. Thus it simplifies the problem to expand the log- 

fa 

arithmic term in 8 in powers of This results in 


R„ = ^ t s _ t o b (1-e 

s * ~ir 


(At s /t 0 ) 


) 


(8b) 
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where only terms up to the first power of b/a have been retained. 

The next step is to solve 6, 7 and 8b simultaneously for a, t 
t 0 . From 8b, 


- (Ato/t 0 ) 

t D b e °' = R s - At s + t 0 b 


( 10 ) 


From 7, 


t c b e " (At s/ t o) - to (t f - a) 

a 2 a 2 


and from 6, 


b = T i - a 


Substitution of these last two expressions in 10 yields: 


R s - + 0° (Ti - T f ) = 0. 

a ^ 


Upon rearrangement this can be written 


i 2 M 
At- 


+ (Ti - T f ) = 0. 

At c 1 r 


From 6 and 7 it is seen that: 


(At s /t 0 ) _ - a 


T i " a 


so that 


or 


In ( Ti ~ a ) = 

Tf - a c o 


At, 


fc o = 


In Cill ) 


T fa 


(11) 


and 


i 
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Substitution of this into 10 yields 


a 


2 Rg 

At s 


T i ' T f 

ln ( T i ~ a ) 
T f - a 


= 0 . 


( 12 ) 


It is possible to find an approximate solution to 12 for a by use 
of Newtons method [see, for example, Franklin pages 26-28], This is an 
iterative method for finding the roots of an equation of the form f(X) = 0. 
In essence, an initial estimate of the root X Q is used to find an improved 
estimate from the equation 

i df(X) 

x i " x ° " laj • 


Thus in the present problem one can write 
f(a) = 0 where f(a) = 


Then 


a 2 R, 


At, 


s ' i 

— - a + 


T « - T f 


in ( T i~ a ) 
T f -a 


df = 2a R s - 1 - 


da 


At, 


~ T i ' T £ 


ln( T i' a )| (Tf-aHTf-a) 

T f-a 


Once a has been determined this way it is a simple matter to use 6 


to determine b and 11 to obtain t Q . 


Now an integration of 1 yields the actual number of rotations of 


the satellite, i.e., 


t iS At ± 


R(t) = 


Rl + t-Atj 


to t> (1 - exp - 


J ), At i ^t<L 1 


t-Li 


R 2 + 


T, + 


1^ £ t £ L 2 


OCt-Lp 


where R^ = At^/r^ 
and 

R2 = R i + (1 - exp { - }) 

1 a t o 

Application 

An application of the method described above to the pass shown in 
Figure 2 is useful for illustrative purposes. The input parameters are 
as follows: 

Tj> = 2.2935 seconds 

Tf = 2.2909 seconds 
At s = 30.5 sequences 

R s = 1088.675 revolutions 

The resulting values for a, b and t c are: 
a = 2.2898 seconds 

b = .003641 seconds 

t Q = 23.96 sequences 

With the additional information that 
t l = 2.29125 

t 2 = 


2.29151 
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and a knowledge of the spin period actually used by the spacecraft, the 
correction angles can be calculated as follows. If 




where t u is the spin period used in the initial analysis, then the 
angle between the apparent sun direction (represented by the psuedo-sun 
pulse) and the true sun direction is given by 

t|r (t)= [R( t ) - R u (t)] rotations 
= 360 [R(t) - R u (t)] degrees 

The correction angles for the pass shown in Figure 2 are listed 
in Table 1* If it is assumed that the magnetic field remained constant 
in direction during this shadow pass, then the apparent direction would 
have rotated through the angle f (t) . It is on this assumption that 
Figure 3 was plotted. The solid curve shows the azimuthal angle of the 
interplanetary magnetic field obtained using the psuedo-sun pulse for 
orientation. The dotted curve shows the relative change i|r(t). 

In conclusion then it may be said that the analysis presented here 
is an effective method for correcting the azimuthal orientation of the 
Explorer 35 spacecraft during its passes behind the moon where it can 
no longer use the sun as a reference. The accuracy of the method is 
estimated to be + 5° to 10°, but can only be checked empirically if it 
is assumed that the magnetic field does not change direction. The largest 
error would occur sometime during the shadow, since at the beginning and 
at the end the orientation has been fitted to the true direction obtained 


by real sun sightings. 


t 

TABLE 

R(t) 

1 

R o (t) 

t (t) 

(Sequences) 

(Revolutions) 

(Revolutions) 

(degrees) 

0 

0.0 

0.0 

0 . 

1 

35.6697 

35.6608 

3. 

2 

71.3395 

71,3216 

6 

3 

107.0114 

106.9823 

10 

4 

142.6856 

142.6431 

15 

5 

178.3618 

178.2982 

23 

6 

214.0401 

213.9532 

31 

7 

249.7203 

249 . 6083 

40 

8 

285.4023 

285.2633 

50 

9 

321.0862 

320.9184 

60 

10 

356.7717 

356.5734 

71 

11 

392.4592 

392.2284 

83 

12 

428.1482 

427.8835 

95 

13 

463.8386 

463.5385 

108 

14 

499.5305 

499.1936 

121 

15 

535.2241 

534.8486 

135 

16 

570.9187 

570.5037 

149 

17 

606.6147 

606.1587 

164 

18 

642.3118 

641.8138 

179 

19 

678.0103 

677.4688 

195 

20 

713.7100 

713.1238 

211 

21 

749.4106 

748.7789 

227 

22 

785.1123 

784.4339 

244 

23 

820.8149 

820.0890 

216 

24 

856.5188 

855.7440 

279 

25 

892.2234 

891.3991 

297 

26 

927.9287 

927.0541 

315 

27 

963.6350 

962.7092 

333 

28 

999.3425 

998.3642 

352 

29 

1035.0503 

1034.0192 

371 

30 

1070.7590 

1069.6743 

391 

31 

1106.4687 

1105.3293 

410 

32 

1142.1787 

1140.9844 

430 

33 

1177.8889 

1176.7318 

417 

34 

1213.5979 

1212.4792 

403 

35 

1249 . 3059 

1248.2266 

389 



Figure Captions 


Figure 1 


Figure 2 


Figure 3 


A plot of T-a versus time (solid line) for the recovery of the 
spacecraft following the first shadow pass for which complete 
data were received. Compared with this is the approximation 
derived from equation 1 (dashed line) . 

A plot of spin period versus time for the same shadow pass as 
in Figure 1. The times and spin periods indicated are discussed 
in the text. The solid part of the curve is derived from measured 
spin periods while the dashed portion in the shadow is the approx- 
imation of equation 1 with appropriate values of the constants 
a, b and t Q substituted. 

Change of apparent azimuth of magnetic field observed (solid 
line) during same shadow pass as Figure 1. Compared is the 
predicted change (dashed line) calculated by assuming that the 
magnetic field did not change its direction. 
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